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ADVANCED  ESM  ANGLE  TRACKER 
VOLUME  I  -  THEORETICAL  FOUNDATIONS 


1.  INTRODUCTION 

This  report  summarizes  the  Naval  Research  Laboratory  (NRL)  implementation  of  tracking  filters  used 
for  tracking  bearing,  elevation,  and  electronic  support  measure  (ESM)  parameters.  It  describes  the  theory 
behind  the  NRL  MATLAB  and  C  code  used  for  ESM  angle  tracking,  developed  to  operate  with  the  NRL 
ESM-Advanced  Technology  Demonstration  (ESM-ATD)  and  later  transitioned  to  Lockheed  Martin 
Corporation  for  integration  into  the  Advanced  Integrated  Electronic  Warfare  System  (AIEWS)  program. 
This  code  was  intended  to  act  as  the  ESM  tracker  for  AIEWS  and  to  provide  feedback  control  for  the 
front-end  gating  of  the  ESM  receiver. 

The  Advanced  ESM  Angle  Tracker  was  developed  at  NRL  and  embedded  within  the  Joint  Maritime 
Command  Information  System  (JMCIS)  as  a  prototype  segment.  This  report  presents  the  theory  used  for 
the  tracker  and  represents  Volume  I  of  a  three-volume  set. 

Volume  II  is  a  user’s  guide  [1]  for  running  the  remote  tracker  via  JMCIS.  It  describes  how  users  set 
up  and  control  the  tracker  from  JMCIS  and  display  real-time  track  information,  and  discusses  JMCIS 
Tactical  Database  Manager  (TDBM)  display  integration.  Volume  III  is  a  developer’s  guide  [2]  that 
documents  the  software  development  of  the  Advanced  ESM  Angle  Tracker.  It  describes  the  tracker’s 
graphical  user  interface,  the  TCP/IP  communications  required  to  interface  with  the  remote  process,  and 
the  implementation  of  the  tracker  algorithms  as  a  remote  process. 

The  filters  and  tracking  algorithm  described  in  this  report  are  unique  in  several  ways.  The  bearing 
filter  is  an  Interactive  Multiple  Model  (IMM)  filter  [3,  4]  that  utilizes  three  Kalman  filters  of  different 
orders  to  track  emitter  bearing  angle.  The  requirements  for  combining  first-,  second-,  and  third-order 
tracks  had  not  been  defined  and  implemented  prior  to  this  work.  Another  unique  feature  of  this  tracker  is 
that  it  includes  all  parameters  (bearing,  elevation,  and  ESM  parameters)  in  the  decision  process  for 
correlation  and  association  of  new  reports  with  existing  tracks.  In  addition,  a  new  technique  developed  for 
determining  plant  noise  in  the  ESM  parameters  filter  estimates  plant  noise  as  a  function  of  ESM 
parameter  agility. 

The  IMM  bearing  filter  uses  three  standard  Kalman  filters  as  a  basis  for  filtering  and  predicting 
emitter  bearing  angle.  Filter  selection  for  the  IMM  is  based  on  the  expected  variation  of  emitter 
trajectories  and  the  imposed  spherical  coordinate  system.  The  derivation  of  the  IMM  is  well  documented 
and  is  not  addressed  here;  only  the  unique  NRL  implementation  aspects  of  the  IMM  are  addressed  here. 
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Each  of  the  three  bearing  Kalman  filters  is  tuned  to  a  different  trajectory  model.  The  first  of  the  three 
Kalman  filter  plant  (hereafter  referred  to  as  trajectory)  models  is  a  simple  first-order  filter  that  is  optimal 
for  tracking  stationary  (land-based)  emitters.  This  filter  is  also  suitable  for  tracking  non-maneuvering, 
radially  inbound  or  outbound  platforms.  The  second  Kalman  filter  is  based  on  filtering  both  bearing  and 
bearing  rate.  The  third  Kalman  filter  tracks  bearing,  bearing  rate,  and  bearing  acceleration  so  as  to 
minimize  filter  lag  for  crossing  emitters  and  other  highly  maneuvering  trajectories.  Platforms  tend  to  fly 
in  straight  line  paths;  crossing  platforms  will  produce  an  infinite  number  of  derivatives  in  a  polar 
coordinate  system,  leading  to  significant  filter  lag.  The  defined  second-  and  third-order  Kalman  filter 
trajectory  models  minimize  these  effects. 

Since  it  is  not  known  a  priori  which  model  is  appropriate  for  any  particular  platform  in  track,  the 
NRL  IMM  filter  estimates  the  probability  of  each  trajectory  model’s  applicability  for  the  particular 
platform  under  observation.  These  trajectory  model  probabilities  are  computed  based  on  a  hidden  Markov 
model  (HMM)  as  well  as  on  outputs  from  each  of  the  three  Kalman  filters  used  to  filter  and  predict 
platform  motion. 

The  elevation  filter  is  a  second-order  Kalman  filter  that  tracks  elevation  and  elevation  rate.  The  ESM 
parameters  filters  are  first-order  Kalman  filters  that  can  track  either  constant  or  agile  parameters  by 
automatically  adjusting  the  filter  plant  noise.  The  ESM  parameters  filters  provide  the  basis  for  the  full 
likelihood  association  used  in  associating  intermittent  bearing/elevation  tracks  in  dense  scenarios. 


2.  IMM  BEARING  FILTER  IMPLEMENTATION 

IMM  bearing  filter  update  consists  of  the  following  steps. 

1 .  IMM  bearing  filter  initialization,  which  consists  of  IMM  and  Kalman  filter  initialization. 

2.  The  IMM  track  update  process,  which  consists  of  state  estimator  mixing  (interaction), 
Kalman  filter  updates,  model  probabilities  update,  and  combining  or  merging  of  the  three 
state  estimates. 


2.1  IMM  Bearing  Filter  Initialization 

The  IMM  bearing  filter  assumes  a  hidden  Markov  model  that  controls  trajectory  model  switching 
with  a  switching  probability  that  is  assumed  known.  This  transition  matrix  is  defined  as 

‘\-\Pii\-  (!) 


where  the  transition  probabilities  are 

Pq  =  Prob{Af/(£)|Mi(£-l)},  ij  =  1,2,3,  (2) 

and  M  .  ( k )  stands  for  “model  j  in  effect  during  the  period  ending  at  time  k 

The  above  expression  should  be  interpreted  as  the  probability  that  the  trajectory  model  for  the  track  at 
time  k  equals  model  j  given  that  at  the  previous  update,  the  trajectory  model  was  i  where  (/,/)  range  from 
1  to  3.  For  example,  in  the  NRL  model,  the  three  elements  in  the  first  row  of  the  matrix  indicate  the 
transition  probability  from  the  bearing-only  model  to  the  bearing-only  model  (the  1,1  element),  to  the 
second-order  model  (the  1,2  element),  and  to  the  third-order  model  (the  1,3  element),  respectively.  The 
second  and  third  rows  should  be  interpreted  in  a  similar  manner. 
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The  NRL  IMM  bearing  filter  initializes  this  Markov  model  transition  matrix  as 


Ai 

(1- P22 )/2 

il~Pn)/2 


(!"A  i)/2  il-Pn)/2 
Pi  2  {l~Pll)/2 

il~P3l)/2  P33 


This  matrix  is  initialized  by  setting 


Ail 

".90" 

Pn 

= 

.95 

_P33  \ 

.90 

(3) 


(4) 


The  transition  probability  pn  should  be  inteipreted  as  the  probability  that  the  first-order  model 

transitions  to  itself  (or  stays  a  first-order  model)  during  the  bearing  track  update  process.  Due  to  the  form 
of  the  matrix  specified  above,  there  is  also  a  small  and  equal  probability  that  the  current  first-order 
platform  model  will  transition  to  the  second-order  or  third-order  models  during  track  update.  The  two 
remaining  transition  probabilities,  p22  and  pi3 ,  should  be  interpreted  in  a  similar  manner.  We  note  that 

the  NRL  model  gives  preference  to  the  second-order  bearing  filter  by  setting  the  return  probability  to  a 
slightly  higher  value.  Thus,  a  track  that  is  second-order  (maintains  bearing  and  bearing  rate)  is  more  likely 
to  stay  second-order  during  the  track  update  process  than  a  first-  or  third-order  track. 

The  initial  probability  distribution  associated  with  the  three  trajectory  models  (hereafter  referred  to 
only  as  models)  is  also  assumed  known.  In  the  NRL  model,  this  initial  distribution  is  given  by 


Pi 


Pi 

Pi  > 
P3_ 


(5) 


where  this  vector  is  initialized  to 


Pi 


.05 

.85 

.10 


(6) 


Again,  the  NRL  model  gives  preference  to  initiating  second-order  tracks.  The  assoc_f  routine 
correlates  current  cluster  descriptors  (CDs)  with  prior  CDs  that  have  been  previously  stored  and  that  have 
not  correlated  with  existing  tracks  (defined  as  first  points).  A  CD  is  equivalent  to  an  ESM  contact 
consisting  of  extracted  information  from  a  set  of  post-processed  pulse  descriptor  words  (PDWs).  If  a 
correlation  exists  in  the  appropriate  parameter  space,  a  new  tentative  bearing  track  is  initiated.  Although 
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the  transition  matrix  is  a  constant,  the  model  probabilities  are  updated  on  each  track  update  as  discussed 
below. 


2.2  IMM  Track  Update  Process 

Computations  associated  with  the  IMM  track  update  process  are  discussed  below. 

2.2.1  State  Estimator  Mixing  (Interaction)  and  State  Prediction 

For  any  initiated  track,  once  an  associated  contact  is  identified,  the  first  step  in  the  IMM  bearing 
update  process  is  mixing  the  previous  filter  estimates  stored  in  the  extended  track  file  (ETF).  This  process 
is  known  as  mixing,  or  interaction,  and  allows  the  IMM  to  achieve  the  accuracy  associated  with  a  second- 
order  generalized  pseudo-Bayesian  (GPB2)  filter  [4]  while  maintaining  three  rather  than  nine  separate 
Kalman  filters. 

The  3x3  mixing  probability  matrix  is  given  by 


P  = 

>\j 


C, 

Pv 

c, 


^•A 

C2 

Pl2 

'  P2 

c2 

P32 

~P3 

C, 


Pn 

-‘Pi 


C3 

P3  3 

~P3 


(7) 


where  c  is  a  3  x  1  normalizing  constant  given  by 


c  = 


=  P 1  p 

ij^j 


(8) 


and  P  is  the  complex  conjugate  transpose  of  P- . 


Using  the  above  mixing  probabilities,  the  mixed  filter  estimates  are  obtained  as 


x 


(9) 


where  xj  (k )  is  the  mixed  estimate  for  model  j,  xt  (k )  is  the  Ith  filtered  estimate  (ij  =  1,2,3)  obtained 
from  the  ETF,  and  pt ;  are  the  mixing  probabilities  defined  above.  The  associated  scalar  covariances  are 
given  by 
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vl  (k)  =  ILPij  C k)\vi  (*0+[*.-  (k)-*7  (*)][*<  (k)-*7  (A0] 


(10) 


It  should  be  noted  that  for  j  =  1,  both  x"'  (k )  and  \\  ( k )  arc  1  x  1  or  scalars.  For  j  =  2,3,  this  is  not 

the  case.  Since,  in  general,  the  filter  estimates  vary  in  size  from  1  x  1  for  the  first-order  filter  to  3  x  1  for 
the  third-order  filter,  mixing  the  ETF  estimates  is  not  straightforward.  Before  computing  the  above 
mixing  equations,  it  is  necessary  to  resize  the  filter  estimates  and  their  associated  covariance  to 
appropriate  sizes  for  mixing.  This  can  be  called  upsizing  and  downsizing.  For  example,  assume  that  the 
form  of  the  ETF  estimate  and  covariance  for  the  third-order  filter  may  be  written  as 


and 


( k ) 

4  (aO 

x3  (&)  = 

(A:) 

= 

4w 

x33)  (A;) 

4  (AO 

v3(U)(*) 

vr(k) 


vr(k) 

A2’2)(k) 

vf’2)(*) 


vf’3)(*) 

v^(k) 

vf’3)(*) 


(11) 


(12) 


where  b-,  (k )  is  the  filtered  bearing  estimate.  Downsizing  is  simply  the  elimination  of  the  appropriate  row 
and  column  before  mixing. 


Downsizing  the  third-order  filter  for  mixing  to  obtain  a  mixed  first-order  estimate  is  accomplished  by 
defining  the  intermediate  downsized  third-order  filter  estimate  as 


V31 


(13) 


and  associated  covariance  as 


viAk)= 


(14) 


Upsizing  in  the  NRL  implementation  implies  the  addition  of  zero  elements  to  expand  the  filter 
estimate  and  associated  covariance  for  mixing.  Consider  the  j  =  3  case  above.  Both  the  first-  and  second- 
order  filters  require  upsizing  before  the  mixing  process.  To  upsize  the  first-order  filter,  an  intermediate 
vector  is  defined  as 
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Xj  (A) 
0 
0 


and  the  associated  3x3  upsized  covariance  is  defined  as 


Vf,  «  = 


Vj  (A')  0  0 

0  0  0 
0  0  0 


Similarly,  the  second-order  filter  is  upsized  to 


%(*)= 


x. 


4”  (A) 

(2)  (A) 


(15) 


(16) 


(17) 


and  the  associated  3x3  upsized  covariance  is  defined  as 


v 


U 

23 


(A) 


vf(A)  v?’2)(A)  0 

vf}(A)  vf'2)(A)  0 

0  0  0 


(18) 


All  other  upsizing  and  downsizing  are  handled  in  a  similar  manner. 

2.2.2  Kalman  Filter  Updates 

The  three  Kalman  filters  should  be  treated  as  subroutines  for  the  IMM  bearing  filter.  Therefore, 
Kalman  filter  updates  may  be  performed  by  passing  in  the  bearing  data  and  the  mixed  estimates  xj  (A) 

and  their  covariance  vj  (A)  in  lieu  of  the  normal  ETF  values.  The  corresponding  output  of  the  Kalman 

filters  for  time  A  are  defined  as  x.(A)  and  v  (A),  respectively.  These  outputs  are  used  for  merging  the 

estimates  using  the  updated  model  probabilities.  These  values  are  also  saved  in  the  ETF  along  with  the 
merged  estimate  and  its  covariance.  In  addition,  each  of  the  three  Kalman  filters  sends  back  the 
innovation  sequence  and  its  covariance,  which  are  used  for  the  update  of  model  probabilities. 
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2.2.3  Update  of  Model  Probabilities 

The  model  probabilities  are  updated  in  two  steps.  First,  the  likelihood  associated  with  each  of  the 
filter  updates  is  computed.  This  likelihood  is  defined  as 

hj{k)  =  p{z(h)\M]{k),Zk-1),  (19) 

where  z(&)  is  the  current  bearing  measurement  from  the  associated  contact  and  7f  1  is  a  vector 
representing  all  the  data  through  time  k  —  1 .  This  may  be  computed  as 


A,(*)= 


expj-ibJSjVd} 

(  r  -ni/2 

jdet  2k Sj  | 


(20) 


where  y  .  is  the  innovation  obtained  during  the  update  of  the  /h  Kalman  filter  and  S  ■  is  its  associated 
covariance.  Since  the  dimensionality  of  S .  is  always  equal  to  the  dimensionality  of  the  measurement 
space,  S  .  is  always  a  scalar  and  the  inverse  indicated  above  never  needs  to  be  computed.  Therefore,  the 
above  may  be  simplified  to 


A,«  = 


(21) 


Based  on  the  likelihoods  defined  above,  the  updated  model  probabilities  may  be  computed  as 


Pj  ( k ) 


where  a  second  normalizing  constant  is  defined  as 


c=EA;W'c/ 


j 


(22) 


(23) 


and  Cj  was  previously  defined. 
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2.2.4  Merging  the  State  Estimates 

The  merged  estimates  and  their  covariances  are  obtained  by  computations  similar  to  that  discussed 
above  for  mixing  (Section  2.2.1).  The  merged  minimum  mean-squared  error  (MMSE)  filter  estimate  is 
given  by 


x(k)  =  YJXl(k)pi(k), 

i 


(24) 


and  its  associated  covariance  by 

v{k)  =  ^p^k-\){vXk)  +  \xXk)-x(k)\xi(k)-x(k)]  J,  (25) 


where  the  individual  estimates  and  their  covariance  must  be  upsized  and  downsized  appropriately  as 
discussed  in  Section  2.2.1.  The  correlation  and  association  process  uses  the  estimates  associated  with  the 
individual  filters. 


2.3  Extended  Track  File  Entries 

The  track  file  database  must  save  the  three  separate  Kalman  filter  estimates,  x  .(&),  which  vary  in 
size  from  1  x  1  to  3  x  1,  and  their  associated  covariances,  Vj(k),  which  vary  in  size  from  1  x  1  to  3  x  3. 
In  addition,  the  merged  estimate  x{k )  and  its  covariance  v(k )  should  be  saved. 


2.4  The  Single-State  Bearing  Filter 

The  following  describes  the  mathematics  for  the  recursive,  single-state  bearing  angle  filter.  This 
report  describes  both  the  initialization  and  the  update  of  the  filter.  The  purpose  of  the  filter  is  to  estimate 
bearing  for  stationary  or  radially  inbound  emitters  and  to  compute  optimal  time-varying  bearing  gates. 

The  bearing  filter  is  a  first-order  Kalman  filter  that  models  only  bearing.  This  is  one  of  the  three 
independent  bearing  filters  used  in  the  IMM  architecture.  The  filter  assumes  that  the  measurement  error 
variance  is  properly  defined  for  the  ESM  bearing  angle  as  a  function  of  measured  frequency  and  signal- 
to-noise  ratio  (SNR). 


2.4.1  Overview 

The  bearing  filter  described  is  a  first-order  Kalman  filter.  The  filter's  general  equation  corresponds  to 
an  alpha  filter  and  is  given  by 


Xj  (A')  =  Oj  (£, k  -l).*)  (k-\ )  +  Kx  (&)  z[k)-Hl  (£)©j  {k,k-  l)xj  (£-1)]  ,  (26) 


where  the  subscript  indicates  the  first-order  filter.  For  the  first-order  bearing  filter  we  define 
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Xi(A)  =  [Ai(A)J, 

(27) 

Oj  (A,  A-l)  =  [1], 

(28) 

z(k)  =  bm(k ), 

(29) 

(*)  =  [!]■ 

(30) 

In  the  above,  X,  (A  —  l)  is  the  lxl  vector  of  filtered  bearing  estimates  valid  at  time  ( A'  —1 )  from  the 
previous  iteration  and  valid  before  the  current  bearing  data  is  received.  Time  [k  —  1)  is  the  last  time  the 
track  was  updated  (time  of  previous  cluster  update),  and  time  k  is  the  current  cluster  measurement  time, 
z  (  A )  —  bm  (A  )  is  the  bearing  angle  measurement  from  the  current  CD.  O,  (A,  A  —  1)  is  a  1  x  1  state  or 

plant  transition  model  and  describes  the  emitter  motion  from  time  segment  (A  —  1)  to  k  .  Hl  (k)  is  a  1  x 

1  measurement  transformation  vector,  and  due  to  the  form  of  Hx ,  K  { (k )  is  a  1  x  1  vector  of  Kalman  gain 
coefficients  that  weight  the  measurement  residuals. 

The  above  general  equation  for  the  filter  update  (Eq.  26)  may  also  be  simplified  by  introducing  the 
predicted  estimates  and  the  predicted  measurement  terms.  In  terms  of  the  predicted  estimates,  the  general 
form  of  the  filter  equation  may  be  written  as 


x,  (A)  =  x,  (A,A-l)  +  Kx  (A)[z(A)-//,  (A)x x  (A,  A- 1)]  , 

where 


xx(k, k  -1)  =  Oj  (A,  A  -  l)X[  (A- 1)  =  bx(k,k- 1), 


(31) 


(32) 


and  bx(k,k-\)  is  defined  as  the  predicted  bearing.  Defining  the  predicted  measurement  as 


(k,k- 1)  =  Hx  (A')  Xj  ( k , k  - 1)  =  bx  (A,  k  - 1)  (33) 

allows  the  final  form  of  the  general  Kalman  filter  equation  to  be  written  in  compact  form  as 

Xj  (A)  =  Xj  (A,  A-l)  +  KX  (A)[z(A)-Zj  (A,  A-l)] .  (34) 

Due  to  the  form  of  Hx  above,  the  form  of  the  Kalman  gain  equation  is  recognized  as  1  x  1.  Taking 
this  into  account  and  expanding  the  above  yields  one  simple  filter  equation  for  bearing: 
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bt C k )  =  fe, (k, k-l)  +  Kl  (k)-  bm  ( k)-bx ( k , k  - 1) 


(35) 


The  sequence  yx{k)  =  z{kj  —  zl[k,k  —  i)  is  known  as  the  innovation  sequence  and  represents  the 
new  information  to  the  filter.  This  term  is  also  passed  up  to  the  IMM  bearing  filter  for  model  probability 
updates.  Finally,  expanding  the  predicted  estimate  defined  as  ®j  [k,k-\)xx  (A'  —  1 )  yields  the  expansion 
for  the  filter  predictions  as 


bx(k,k-\)  =  b1  (A-l) .  (36) 

2.4.2  Initialization 

The  single-state  bearing  angle  Kalman  filter  may  be  initialized  using  one  of  two  equivalent 
methods. 

Method  1 

In  this  method,  the  bearing  filter  is  initialized  with  a  large  variance: 


Xj(0)  =  0  and  Vj(o)  =  1012,  (37) 

where  v1  (0)  is  the  initial  lxl  filtered  covariance  value.  A  large  covariance  value  simply  indicates  no  a 
priori  information  and,  therefore,  large  uncertainty. 

Method  2 

In  this  method,  the  bearing  filter  is  initialized  from  the  first  cluster  measurement  of  bearing  angle.  In 
this  case  we  define 


^i(1)  =  Z,m(1)and  vi(1)  =  v™c(1)> 


(38) 


where  bm(  1)  is  the  bearing  measurement  provided  with  the  first  cluster  and  vmc  ( 1 )  is  the  computed 

measurement  variance  associated  with  the  same  bearing  measurement  discussed  in  Section  2.4.5.  The  first 
measurement  is  used  to  initialize  the  filter  and  the  filter  iteration  starts  at  time  k  =  2  with  the  second 
measurement,  while  in  Method  1 ,  the  iteration  starts  at  time  k  =  1  with  the  first  measurement. 


2. 4. 3  Predicted  Estimate 

When  a  new  measurement  is  associated  with  the  bearing  track,  the  first  step  in  the  filter  update 
process  requires  that  the  ETF  estimate  be  predicted  to  the  time  of  the  current  associated  measurement 
(time  alignment).  For  the  bearing  filter,  the  predicted  estimate  is  obtained  from  the  transition  matrix 
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defined  above  (Eq.  28),  the  filtered  estimate  x,  (k  -  \)  from  the  previous  ETF  update,  and  the  time 
difference  between  the  associated  contact  and  the  ETF  time  (time  of  the  last  track  update)  as 


x,  (k,k- 1)  =  Oj  [k,k-  l)xj  (&-l) 


(39) 


or,  equivalently, 


bl{k,k-\)  =  bl(k-l) ,  (40) 

where  the  lxl  vector  Xj  {k,  k  —  1 )  is  the  prediction  of  x,  (k  —  1 )  to  the  current  CD  time  and  xt  (t-i)  is 
the  previous  filtered  estimate  obtained  from  the  ETF. 


2. 4. 4  Predicted  Variance  and  Plant  Noise  Estimates 

Once  the  prediction  is  computed,  computation  of  the  Kalman  filter  gain  is  required  to  merge  the 
bearing  CD  data  with  the  bearing  angle  of  a  track.  The  Kalman  gain  computation  requires  an  estimate  of 
the  predicted  variance  and  the  measurement  variance.  The  lxl  predicted  variance  is  given  by 


V;  (k ,  k  - 1)  =  Vj  (k  - 1)  +  Vj  ( k ) , 


(41) 


where  v)  ( k )  is  the  lxl  plant  noise  estimate  and  Vj  (k  —  l)  is  obtained  from  the  track  file  (ETF).  For  the 
bearing  filter,  the  plant  noise  is  set  to 


where 

b[a)  (&)  =  1. 


(42) 


(43) 


w 


is  a  constant  and  should  be  approximately  set  to 


w\a)  =0.1;r/ 180/ 40  ^0.00005. 


(44) 


Therefore, 


v,  (k)  =  0.00005 . 


(45) 
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As  a  final  note,  w\a^  is  a  tuning  parameter  that  controls  filter  bandwidth  and  should  be  optimized  during 
testing. 


2.4.5  Filter  Gain  and  Measurement  Noise  Estimate 
The  lxl  Kalman  gain  Kx  is  given  by 


K,(k) 


v,  (k,k-\)H[ 
Hlvl(k,k~l)H[  +  vmc(k)' 


(46) 


We  note  here  that  a  simple  divide  is  used  and  that  no  matrix  inverse  operation  is  required  since  the 
term  Hlvl  {k,k  —  \)H[  is  a  scalar.  However,  to  compute  this  gain,  an  estimate  of  CD  bearing 

measurement  noise  is  required.  If  the  CD  bearing  measurement  noise  is  not  included  in  the  cluster 
message  it  can  be  computed  from  the  single  pulse  specification  variance  according  to  the  following 
equation: 


vmc(k)  = 


V,n  (k) 


(47) 


where  the  variable  vmc  (&)  is  the  measurement  noise  variance  in  bearing  associated  with  the  cluster 
estimate.  The  variable  vm  (k  )  is  the  single  pulse  measurement  accuracy  (equivalent  to  the  PDW 

accuracy)  of  the  bearing  angle,  and  n  is  the  number  of  pulses  reported  in  the  CD  message.  For  the  bearing 
estimate,  the  measurement  accuracy  is  a  function  of  frequency,  SNR,  and  measured  emitter  bearing  angle. 
These  features  must  be  represented  as  such  in  a  suitable  tabular  format.  This  representation  should  be 
obtained  from  system  testing  and  should  follow  the  form 


■W  = 


/o 


/»(*) 


(48) 


where  <7  a  should  be  characterized  and  is  the  base  bearing  accuracy,  which  is  a  function  of  measured 

emitter  bearing  angle,  and  f0  is  the  constant  frequency  at  which  <Ja  is  measured.  The  term  fm  ( k )  is  the 

frequency  measurement  associated  with  the  bearing  angle  provided  in  the  cluster  descriptor.  vm  (k  ) 

should  be  inversely  proportional  to  SNR  and  limited  to  a  constant  that  represents  best-case  bearing 
accuracy. 

The  actual  value  of  <7a  must  be  obtained  from  system  measurements  and  testing  and  used  here  if 

reliable  tracking  is  to  be  expected.  This  same  comment  applies  to  the  elevation  filters.  Note  also  that  cra 

should  be  an  over-bounded  approximate  estimate  of  accuracy  in  bearing  angle  as  a  function  of  frequency 
and  measured  bearing  angle. 
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The  Kalman  gain  may  also  be  written  in  terms  of  the  covariance  of  the  innovation  sequence  as 


Kt(k) 


v,  (4,4-1  )H[ 


(49) 


where  S\(k)  is  the  scalar  covariance  of  the  innovation  sequence  defined  above  and  is  given  by  the 
equation 

Si  (k)  =  //,v,  ( k ,  k  - 1  )H[  +  vmc  (A) .  (50) 


This  term  is  also  passed  up  to  the  IMM  bearing  filter  for  model  probability  updates. 


2. 4. 6  The  Filtered  Estimate 

After  the  Kalman  gains  are  computed,  the  new  filtered  or  ETF  track  estimate  is  given  by 


Xj  (A)  =  xx  (k,k-\)  +  Kx  (A){z(A)-Zj  (A,A-l)j  (51) 

or 

x,  (k)  =  x,  ( k , k-l)  +  Kl  (&) yx  (k) ,  (52) 

where  xx  (A)  is  the  filtered  bearing  estimate,  z (A)  is  the  current  bearing  cluster  measurement,  and 
yl  (A)  is  the  innovation  sequence.  As  mentioned  above  (Eq.  35),  this  is  equivalent  to  computing 


bx  (A)  =  bt  (A,  k-l)  +  Kl  (A)  yx  (A)  . 


(53) 


2.4.7  Filtered  Variance 

The  filtered  variance  of  this  estimate  is  simply 


v1(A)  =  [l-K1(A)//i]vi(A,A-l), 


(54) 


which  may  be  expanded  as 


Vj  (A)  =  l-Kj(A)  Vj(A,A-l). 


(55) 
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2. 4. 8  Track  Coasting 

There  will  be  times  when  an  emitter  bearing  estimate  is  not  available.  For  example,  this  may  be  due 
to  ESM  system  failures,  improper  dwell  scheduling,  and  emitter  fading  due  to  low  SNR.  Since,  in  these 
cases,  bearing  estimates  will  not  be  available,  the  bearing  track  should  be  coasted  to  the  time  associated 
with  the  bearing  update.  The  coasting  duration  time  is  a  function  of  track  confidence  as  well  as  track 
status.  This  is  discussed  in  a  different  document. 

If  the  bearing  ETF  track  is  in  coast  mode,  the  corresponding  elevation  and  ESM  parameters  must  also 
be  in  coast  mode.  It  is  possible,  and  at  times  necessary,  to  update  a  bearing  track  while  coasting  the 
corresponding  elevation  track. 

It  should  also  be  noted  that  for  these  reasons,  the  initiate  and  drop  track  criteria  for  bearing  angle 
should  be  evaluated  separately  from  the  associated  elevation  track.  That  is,  bearing  tracks  should  be 
allowed  to  initiate  and  drop  without  reference  to  the  elevation  track.  However,  as  mentioned  earlier,  if  a 
bearing  track  does  not  continue  to  exist,  then  the  corresponding  elevation  and  ESM  parameter  track  will 
also  cease  to  exist. 

The  simplest  way  to  derive  the  filter  for  coasting  is  to  recognize  that  coasting  is  equivalent  to 
receiving  a  measurement  with  zero  accuracy  (or  infinite  measurement  variance).  Using  this  analogy 
implies  the  Kalman  gain  (see  the  appropriate  equation  above)  goes  to  zero  and  the  filtered  estimate  simply 
becomes  the  predicted  estimate.  The  sequence  for  filtering  during  a  bearing  coast  is  the  same  for  each  of 
the  three  bearing  filters.  Coasting  a  track  is  achieved  by  setting  the  filter  gain  to  zero  (equivalent  to  setting 
the  measurement  variance  to  infinity). 

With  the  Kalman  gain  set  to  zero,  the  remaining  equations  for  filtered  estimate  and  variance  are  the 
same  as  those  defined  above  (Eqs.  51-53).  In  this  case,  the  filtered  estimate  is  given  by  the  prediction 


jCj  (k)  =  x,  (k,  k  - 1) 


(56) 


or,  equivalently, 


b](k)=bt(k,k-\),  (57) 

and  the  filtered  variance  takes  the  form 

vl(k)=vl(k,k-l).  (58) 


2.5  The  Two-State  Bearing  Filter 

The  following  describes  the  theory  used  for  the  recursive,  two-state  bearing  angle  filter.  This 
section  describes  both  the  initialization  and  the  update  of  the  filter.  The  purpose  of  the  filter  is  to  estimate 
bearing  for  nominal  trajectory  emitter  platforms  and  to  compute  optimal  time- varying  bearing  gates. 


Advanced  ESM  Angle  Tracker,  Volume  I 


15 


The  bearing  filter  is  a  second-order  Kalman  filter  that  models  both  bearing  and  bearing  rate.  This 
is  one  of  the  three  independent  bearing  filters  used  in  the  IMM  architecture.  The  first-  and  third-order 
bearing  angle  Kalman  filters  are  described  in  separate  sections.  The  filter  assumes  that  the  measurement 
error  variance  is  properly  defined  for  the  ESM  bearing  angle  as  a  function  of  measured  frequency  and 
SNR. 


2.5.1  Overview 

The  bearing  filter  described  is  a  second-order  Kalman  filter.  The  filter's  general  equation  corresponds 
to  an  alpha-beta  filter  and  is  given  by 

x2  (&)  =  ®2  (k,k  —  l)x2  {k-\)  +  K2  (£)[z(£)-//2  (£)®2  (k,k  —  l)x2  (£-1)]  ,  (59) 

where  for  the  bearing  filter  we  define 


x2  (k) 


b2(k) 

k(k) 


(60) 


and 


®2{k,k- 1)=  Q 

z(k)  =  bm(k), 

H2(k)  =  [\  0], 


(61) 


(62) 


(63) 


In  the  above,  x2  (k  —  l)  is  the  2x1  vector  of  filtered  bearing  estimates  valid  at  time  k  —  I  from  the 

previous  iteration  and  valid  before  the  current  bearing  data  is  received.  At  is  the  difference  between  the 
current  cluster  measurement  time  and  the  last  time  the  track  was  updated  (time  of  previous  cluster 
update).  z{k}  =bm  (k  )  is  the  bearing  angle  measurement  from  the  current  CD.  ®2  (k,k-\)  is  a  2  x  2 

state  or  plant  transition  model  and  describes  the  emitter  motion  from  time  segment  k  —  1  to  k  .  II 2  (k )  is 
a  1  x  2  measurement  transformation  vector,  and  due  to  the  form  of  H2  ,  K2{k )  is  a  2  x  1  vector  of 
Kalman  gain  coefficients  that  weight  the  measurement  residuals. 

The  above  general  equation  for  the  filter  update  may  also  be  simplified  by  introducing  the  predicted 
estimates  and  the  predicted  measurement  terms.  In  terms  of  the  predicted  estimates,  the  general  form  of 
the  filter  equation  may  be  written  as 


x2  (£)  =  x2  {k,k  -\)  +  K2  (k)  z(k)-H2  (A:)x2  (&,&-!)] 


(64) 
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where 


x2(k,k- 1)  =  ®2  {k,k- l)x2  (k- 1) 


^(A:,A:-1) 

b2  [k-l)  +  At-b2  (k- 1) 

b0  (k,  k  - 1) 

b2(k-l) 

(65) 


where  b2  ( k ,  k  - 1)  is  defined  as  the  predicted  bearing  estimate  and  b2  (k,  k  —  1)  is  defined  as  the  predicted 
bearing  rate  estimate.  Defining  the  predicted  measurement  as 


z2  [k,k- 1)  =  H2  (£)x2  [k,k- 1)  =  b2  (k,k  -1)  (66) 

allows  the  final  form  of  the  general  Kalman  filter  equation  to  be  written  in  compact  form  as 

x2  (£)  =  x2  (k,k-\)  +  K2  (&)  jz(&)-z2  (k,k- l)j  .  (67) 

Due  to  the  form  of  H2  above,  the  form  of  the  Kalman  gain  is  recognized  as  2  x  1.  Taking  this  into 

account  and  expanding  the  above  yields  two  simple  filter  equations  for  the  bearing  and  bearing  rate 
filters: 


b2(k)  =  b2(k,k  -\)  +  K2^  ( k)\z(k)-z2  (k,k- 1)} 
b2 (k)  =  b2(k,k  -l)  +  (k){z(k)~ z2  [k,k- 1)}. 


(68) 


The  sequence  y2(k}  =  z{k}  —  z^ik.k  —  is  known  as  the  innovation  sequence  and  represents  the 

new  information  to  the  filter.  This  term  is  also  passed  up  to  the  IMM  bearing  filter  for  model  probability 
updates. 

Finally,  expanding  the  predicted  estimate  defined  as  ®2  (k,k-  l)x,  (&- 1)  yields  the  expansion  for 
the  filter  predictions  as 


b2{k,k- 1) 

b2  [k  -  \)  +  At  -b2  (k  —  l) 

b2(k,k-\) 

b2(k- 1) 

(69) 


2.5.2  Initialization 

The  two-state  bearing  angle  Kalman  filter  may  be  initialized  using  any  one  of  three  equivalent 
methods. 
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Method  1 

In  this  method,  the  bearing  filter  covariance  matrix  is  initialized  with  large  values  on  the  diagonal: 

v2(0)  =  1012I2,  (70) 

0" 

1 

v2  (0)  is  the  initial  2x2  filtered  covariance  value  and  I2  is  a  2  x  2  identity  matrix.  A  large  value  simply 
indicates  no  a  priori  information  and,  therefore,  large  uncertainty. 

Method  2 

In  this  method,  the  bearing  filter  is  initialized  from  the  first  cluster  measurement  of  bearing  angle.  In 
this  case  we  define 


where 


x2(0)  = 


and 


I2  = 


UD 

o 


and  v2  (1)  = 


Vmc  (0  0 

a  i  nl2 


(71) 


where  bm(  1)  is  the  bearing  measurement  provided  with  the  first  cluster  and  vmc  ( 1 )  is  the  computed 

measurement  variance  associated  with  the  same  bearing  measurement.  The  filter  iteration  starts  at  time 
k  =  2  with  the  second  measurement. 

Method  3 

In  this  method,  the  bearing  filter  is  initialized  from  the  first  two  cluster  measurements  of  bearing 
angle.  In  this  case,  we  define 


and 


b,n(  2) 

(2) -£,„(!) 

At 


vi(2) 


Vmc  (2) 
Vmc  (2) 

At 


Vn,c  (2) 

At 

2v  (2) 


(72) 


(73) 


The  filter  iteration  starts  at  time  k  =  3  with  the  third  measurement. 
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2. 5. 3  Predicted  Estimate 

When  a  new  CD  is  associated  with  the  bearing  track,  the  first  step  in  the  filter  update  process  requires 
that  the  ETF  estimate  be  predicted  to  the  time  of  the  current  associated  CD  (time  alignment).  For  the 
bearing  filter,  the  predicted  estimate  is  obtained  from  the  transition  matrix  defined  above,  the  filtered 
estimate  x2{k  - 1)  from  the  previous  ETF  update,  and  the  time  difference  between  the  associated  CD  and 
the  ETF  time  (time  of  the  last  track  update)  as 


x2  (k,k- 1)  =  ®2  (k,k-  \)x2  [k- 1) 


(74) 


or,  equivalently, 


b2{k,k- 1) 

b2  (&-l)  +  At-b2  (£-1) 

b2{k,k- 1) 

b2(k- 1) 

(75) 


where  the  2x1  vector  x2  ( k ,  k  —  l)  is  the  prediction  of  x2  ( k  -  l)  to  the  current  CD  time  and  x1  (k  —  l)  is 
the  previous  filtered  estimate  obtained  from  the  ETF. 


2. 5. 4  Predicted  Variance  and  Plant  Noise  Estimates 

Once  the  prediction  is  computed,  computation  of  the  the  Kalman  filter  gain  is  required  to  merge  the 
bearing  CD  data  with  the  bearing  angle  of  a  track.  The  Kalman  gain  computation  requires  an  estimate  of 
the  predicted  variance  and  the  measurement  variance.  The  2x2  predicted  variance  is  given  by 

v2  {k,k-\}  =  v2  (&-l)  +  v2  (&),  (76) 

where  v2  is  the  2x2  plant  noise  estimate  and  v2  ( A:  —  I )  is  obtained  from  the  track  file  (ETF).  For 
the  bearing  filter,  the  plant  noise  is  set  to 

v2(*)  =  bia)  (k)-  w (2a)  ■  b{2]  (k)\  (77) 


where 

At 

1 


(78) 


is  a  constant  and  should  be  approximately  set  to 


w{a)  =0.  Dr/ 180/ 40  ^0.00005. 


(79) 
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Therefore, 


(A:)  =  0.00005  • 


At2 

At 


At 

1 


(80) 


As  a  final  note,  vJf*  is  a  tuning  parameter  that  controls  filter  bandwidth  and  should  be  optimized  during 
testing. 


2.5.5  Filter  Gain  and  Measurement  Noise  Estimate 
The  2x1  Kalman  gain  K2  is  given  by 


K,(k) 


v2(k,k-\)H’2 

H2v2(k,k-\)H'1+vmc(k)' 


(81) 


We  note  here  that  a  simple  divide  is  used  and  that  no  matrix  inverse  operation  is  required  since  the 
term  H2v2  {k,k  -  l)H'2  is  a  scalar.  However,  to  compute  this  gain,  an  estimate  of  CD  bearing 

measurement  noise  is  required.  If  the  CD  bearing  measurement  noise  is  not  included  in  the  cluster 
message  it  can  be  computed  from  the  single  pulse  specification  variance  according  to  the  following 
equation: 


=  (82) 

n 

where  the  variable  vmc(A:)  is  the  measurement  noise  variance  in  bearing  associated  with  the  cluster 

estimate.  The  variable  vm  (k)  is  the  single  pulse  measurement  accuracy  (equivalent  to  the  PDW 
accuracy)  of  the  bearing  angle,  and  n  is  the  number  of  pulses  reported  in  the  CD  message.  The 
computation  of  vm{k)  is  discussed  in  Section 2.4.5. 

The  Kalman  gain  may  also  be  written  in  terms  of  the  covariance  of  the  innovation  sequence  as 


K2  (1) 


v2(k,k-\)H[ 

S2(k) 


(83) 


where  S2  ( k )  is  the  scalar  covariance  of  the  innovation  sequence  defined  above  and  is  given  by  the 
equation 


S2(k)  =  H2v2(k,k-\)H'2+vmc(k) . 


(84) 
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2. 5. 6  The  Filtered  Estimate 

After  the  Kalman  gains  are  computed,  the  new  filtered  or  ETF  track  estimate  is  given  by 

x2  (A')  =  x2  (k, k-\)  +  K2  (&){z(&)-z2  (k,k  -l)| ,  (85) 

where  x2  is  the  filtered  bearing  estimate  and  z(k)  is  the  current  bearing  cluster  measurement.  As 
mentioned  above,  this  is  equivalent  to  computing 


b2(k)  =b2(k,k-l)  +  (k){z(k)-z2  (k,k- 1)} 

b2(k)  =  b2(k,  k  -1)  +  (£){z(A:)-z2  (£,£-!)} 


(86) 


or,  equivalently, 


b2(k)  =  b2(k,k-\)  +  K2^  (k^y2  ( k ) 
<J  ^ 

b2 ( k )  =  b2  ( k ,  k  - 1)  +  K2^  (&)  y2  , 


where  y2  (k  j  is  the  innovation  sequence. 


2.5. 7  Filtered  Variance 

The  filtered  variance  of  this  estimate  is  simply 


v2  (^)  =  [A  ~K2  (A:)//2]  ■  v2  (k,k- 1) ,  (88) 

which  may  be  expanded  as 


v2(k) 


1  -  K2  ]  (k) 
-K?\k) 


■v2{k,k- 1), 


(89) 


where  /2  is  the  2x2  identity  matrix,  and  v2  {k,k  —  1)  was  defined  above. 


2. 5. 8  Track  Coasting 


As  discussed  in  Section  2.5.5,  the  filter  gain  is  set  to 
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K,(k)  = 


0 

0 


(90) 


With  the  Kalman  gain  set  to  zero,  the  remaining  equations  for  filtered  estimate  and  variance  are  the 
same  as  those  defined  above.  In  this  case,  the  filtered  estimate  is  given  by  the  prediction 


x2(A)  =  x2(k,k  - 1) 


(91) 


or,  equivalently, 


b2  ( k )  =  b2  ( k ,  k  —  1) 

(92) 

b2{k )  =  b2(k,k  —  1), 

and  the  filtered  variance  takes  the  form 

v2(k)  =  v2(k,k- 1).  (93) 


2.6  The  Three-State  Bearing  Filter 

The  following  describes  the  theory  for  the  recursive,  three-state  bearing  angle  filter.  This  section 
describes  both  the  initialization  and  the  update  of  the  filter.  The  purpose  of  the  filter  is  to  estimate  bearing 
for  highly  maneuvering  emitter  platforms  and  to  compute  optimal  time-varying  bearing  gates. 

The  bearing  filter  is  a  third-order  Kalman  filter  that  models  bearing,  bearing  rate  and  bearing 
acceleration.  This  is  one  of  the  three  independent  bearing  filters  used  in  the  IMM  architecture.  The  first- 
and  second-order  bearing  angle  Kalman  filters  are  described  in  separate  sections.  The  filter  assumes  that 
the  measurement  error  variance  is  properly  defined  for  the  ESM  bearing  angle  as  a  function  of  measured 
frequency  and  SNR. 


2.6.1  Overview 

The  bearing  filter  described  is  a  third-order  Kalman  filter.  The  filter's  general  equation  corresponds  to 
an  alpha-beta-gamma  and  is  given  by 


x3  (A)  =  0,  (A,  k  -l)x3  (A-1)  +  /l3  (A){z(A)-//3  (A)03  (A,  A  -l)x3  (A  -l)j  ,  (94) 


where  for  the  bearing  filter  we  define 
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4  (*) 

X3(*)=  k(k)  ,  (95) 

_£(*)_ 

"l  At  At2  /l 

®3(k,k- 1)=  0  1  A?  ,  (96) 

0  0  1 

z(k)  =  bm(k),  (97) 

and 

iJ3(*)  =  [l  0  0],  (98) 

In  the  above,  x3  (&  - 1)  is  the  3x1  vector  of  filtered  bearing  estimates  valid  at  time  k  —  1  from  the 

previous  iteration  and  valid  before  the  current  bearing  data  is  received.  At  is  the  difference  between  the 
current  cluster  measurement  time  and  the  last  time  the  track  was  updated  (time  of  previous  cluster 
update).  z(k')  =  bm(k')  is  the  bearing  angle  measurement  from  the  current  cluster  descriptor  (CD). 

0,  (k,  k  -  l)  is  a  3  x  3  state  or  plant  transition  model  and  describes  the  emitter  motion  from  time  segment 
(&  —  1 )  to  k  .  If  3  (k  )  is  a  1  x  3  measurement  transformation  vector  and  due  to  the  form  of  H3 ,  K3  (k)  is 
a  3  x  1  vector  of  Kalman  gain  coefficients  that  weight  the  measurement  residuals. 

The  above  general  equation  for  the  filter  update  may  also  be  simplified  by  introducing  the  predicted 
estimates  and  the  predicted  measurement  terms.  In  terms  of  the  predicted  estimates,  the  general  form  of 
the  filter  equation  may  be  written  as 

x3  ( k )  =  x3  (k,k-\)  +  K3  (A)|z(A)-//3  (£)x3  (k,k- l)j ,  (99) 

where 

b3(k,k-l) 

x3(k,k- 1)  =  ®3  (£,  k  -l)x3  (k- 1)  =  b3(k,k-V) 

b3(k,k- 1) 

4(£-l)  +  Ac4(£-l)  +  0.5-  At2-S3(k-l) 

=  b3  (A:-l)+  At-b3  {k -l) 

k(k- 1) 


(100) 
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where  b, (k, k  - 1)  is  defined  as  the  predicted  bearing  estimate,  b3  ( k , k  —  1)  is  defined  as  the  predicted 

bearing  rate  estimate,  and  b3(k,k  —  Y)  is  defined  as  the  predicted  bearing  acceleration  estimate.  Defining 
the  predicted  measurement  as 


z3  (k,k  —  1)  =  H3(k^x3{k,k-\)  =  b3(k,k-\)  (101) 

allows  the  final  form  of  the  general  Kalman  filter  equation  to  be  written  in  compact  form  as 

x3  (A:)  =  x3  ( k,k -l)  +  K3  (k}{z(k)-z3  [k,k- 1)} .  (102) 

Due  to  the  form  of  the  H3  matrix,  the  Kalman  gain  is  recognized  as  a  3  x  1  matrix.  Taking  this  into 

account  and  expanding  the  above  yields  three  simple  filter  equations  for  the  bearing,  bearing  rate  and 
bearing  acceleration  filters: 


b3(k)  =  b3(k,k-l)  +  K3^  (A:){z(A:)-z3  (k,k  -l)} 

<  b3(k)  =  b3(k, k - 1)  +  K\2^  (k)  jz (A:) - z3  ( k, k - 1)}  (103) 

b3(k)  =  b3(k,k-V)  +  K^  (k){z(k)-z3  (k,k- 1)}. 


The  sequence  y3(k)  =  z(k^)  — z3(k,k  - 1)  is  known  as  the  innovation  sequence  and  represents  the 

new  information  to  the  filter.  This  term  is  also  passed  up  to  the  IMM  bearing  filter  for  model  probability 
updates. 

Finally,  expanding  the  predicted  estimate  defined  as  0(A:,A:-l)x(A:-l)  yields  the  expansion  for 
the  filter  predictions  as 


b3(k,k- 1) 

b3  (k  - 1)  +  A*  •  4  (k  - 1)  +  0.5  •  At2  ■  S3  (k  - 1) 

b3{k,k-\) 

= 

b3  (k-\)  +  At -b3  (A:  - 1) 

b3{k,k- 1) 

L(k-\) 

(104) 


2.6.2  Initialization 


The  three-state  bearing  angle  Kalman  filter  may  be  initialized  using  one  of  three  methods. 
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Method  1 

In  this  method,  the  bearing  filter  covariance  matrix  is  initialized  with  large  values  on  the  diagonal: 


0 

0 

0 


and  v3  (0)  =  1012  I3, 


(105) 


where 


I3 


1  0  0 
0  1  0 
0  0  1 


v3  (0)  is  the  initial  3x3  filtered  covariance  value  and  I3  is  a  3  x  3  identity  matrix.  A  large  value  simply 
indicates  no  a  priori  information  and,  therefore,  large  uncertainty. 

Method  2 

In  this  method,  the  bearing  filter  is  initialized  from  the  first  cluster  measurement  of  bearing  angle.  In 
this  case,  we  define 


w 

'v(i)  0  0  ' 

x3 (1)  = 

0 

and  v3  (l)  = 

0  1012  0 

0 

0 

0 

0 

to 

(106) 


where  bm(  1)  is  the  bearing  measurement  provided  with  the  first  cluster  and  vmc  ( 1 )  is  the  computed 
measurement  variance  associated  with  the  same  bearing  measurement.  The  filter  iteration  starts  at  time 

k  =  2. 

Method  3 

In  this  method,  the  bearing  filter  is  initialized  from  the  first  two  cluster  measurements  of  bearing 
angle.  In  this  case,  we  define 


b„X  2) 

bm(2)-bmQ) 


At 

0 


(107) 
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and 


v 


3 


(2) 


Vmc  (2) 
Vmc  (2) 

At 

0 


U2)  Q 

At 

2v  (2) 

mc  V  o 

(AO2 

0  1012 


(108) 


The  filter  iteration  starts  at  time  k  =  3  with  the  third  measurement. 


2. 6. 3  Predicted  Estimate 

When  a  new  CD  is  associated  with  the  bearing  track,  the  first  step  in  the  filter  update  process  requires 
that  the  ETF  estimate  be  predicted  to  the  time  of  the  current  associated  CD  (time  alignment).  For  the 
bearing  filter,  the  predicted  estimate  is  obtained  from  the  transition  matrix  defined  above,  the  filtered 
estimate  x3  (k  —  l)  from  the  last  ETF  update,  and  the  time  difference  between  the  associated  CD  and  the 
ETF  time  (time  of  the  last  track  update)  as 


x3  {k,k- 1)  =  ®3  {k,k- l)x3  (£-1) 


(109) 


or,  equivalently, 


b^{k,k- 1) 

4(£-l)  +  Ac4(£-l)  +  0.5-At2-4(£~l) 

b^(k,k-\) 

= 

4  [k- 1)  +  AC  4  (k- 1) 

b^(k,k- 1) 

L(k-\) 

(HO) 


where  the  3  x  1  vector  x3  {k,  k  - 1)  is  the  prediction  of  x3  (k  - 1)  to  the  current  CD  time  and  x3  (k  - 1)  is 
the  previous  filtered  estimate  obtained  from  the  ETF. 


2. 6. 4  Predicted  Variance  and  Plant  Noise  Estimates 

Once  the  prediction  is  computed,  computation  of  the  the  Kalman  filter  gain  is  required  to  merge  the 
bearing  CD  data  with  the  bearing  angle  of  a  track.  The  Kalman  gain  computation  requires  an  estimate  of 
the  predicted  variance  and  the  measurement  variance.  The  3x3  predicted  variance  is  given  by 
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v3(k,k-\)  =  v3(k-\)  +  v3(k),  (111) 

where  v3  (A:)  is  the  3x3  plant  noise  estimate  and  v3  {k  - 1)  is  obtained  from  the  track  file  (ETF).  For  the 
bearing  filter,  the  plant  noise  is  set  to 


where 


v3  ( k )  =  b^  (k)-  w3‘-  ■  b3'^  ( k ) 


(fl)  .  d“)  i 


At1/! 

At 

1 


w: 


ia) 


is  a  constant  and  should  be  approximately  set  to 


(112) 


(113) 


w{3a)  =  0.  Dr  7180/40  =  0.00005. 


(114) 


Therefore, 


v3  (A:)  =  0.00005- 


At4/ 4 
At4/! 
At2/! 


At4/! 

At2 

At 


At2/! 

At 

1 


(115) 


As  a  final  note,  is  a  tuning  parameter  that  controls  filter  bandwidth  and  should  be  optimized  during 
testing. 


2. 6. 5  Filter  Gain  and  Measurement  Noise  Estimate 
The  3x1  Kalman  gain  K3  is  given  by 


Ks(k) 


H^(k,k-\)H[+v,c(ky 


(116) 


Note  here  that  a  simple  divide  is  used  and  that  no  matrix  inverse  operation  is  required  since  the  term 
H3v3  (k,k-\)H3  is  a  scalar.  However,  to  compute  this  gain,  an  estimate  of  CD  bearing  measurement 

noise  is  required.  If  the  CD  bearing  measurement  noise  is  not  included  in  the  cluster  message  it  can  be 
computed  from  the  single  pulse  specification  variance  according  to  the  equation 
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:(*)=■ 


w 


(117) 


where  the  variable  vmc(A)  is  the  measurement  noise  variance  in  bearing  associated  with  the  cluster 
estimate.  The  variable  vm  (A )  is  the  single  pulse  measurement  accuracy  (equivalent  to  the  PDW 
accuracy)  of  the  bearing  angle,  and  n  is  the  number  of  pulses  reported  in  the  cluster  descriptor  message. 
The  computation  of  vm(k )  is  discussed  in  Section 2.4.5. 

The  3x1  Kalman  gain  may  also  be  written  in  terms  of  the  covariance  of  the  innovation  sequence  as 


K,(k) 


v3(k,k-l)H'3 

S&) 


(118) 


where  S3  (Aj  is  the  scalar  covariance  of  the  innovation  sequence  defined  above  and  is  given  by  the 
equation 

S3  (k)  =  H3v3  ( k , k  -\)H'3+  vmc  (k) .  (119) 


2. 6. 6  The  Filtered  Estimate 

After  the  Kalman  gains  are  computed,  the  new  filtered  or  ETF  track  estimate  is  given  by 

x,  (A')  =  x3  (A,A-l)  +  K3  (A)|z(A)-z3  (A,  A-l)| ,  (120) 

where  x3(A)  is  the  filtered  bearing  estimate  and  z  ( A )  is  the  current  bearing  cluster  measurement.  As 
mentioned  above,  this  is  equivalent  to  computing 

4(A)  =  4(A,A-1)  +  Af  (A)|z(A)  -z3  (A,  A  - 1)} 

<  4(A)  =  4(A,A-1)  +  Kf  (A)|z(A)-z3 (A, A  — 1)|  (121) 

4  (A)  =  4  (k,  A  - 1)  +  K.f  ]  (A)|z(A)-z3  (A,  A  — 1)| 


or,  equivalently, 
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b3(k)  =  b3(k,k-l)  +  Kf]  (k)y3(k) 

<i3(k)  =  b3(k,k-\)  +  Kf]  (k)y3(k)  (122) 

b3  (k)  =  b3  (k,k  - 1)  +  K®  (k)y3(k), 

where  yi(k)  is  the  innovation  sequence. 


2. 6. 7  Filtered  Variance 

The  filtered  variance  of  this  estimate  is  simply 


v3(k)=  I3-K3{k)H3  v3{k,k- 1), 


which  may  be  expanded  as 


v 


3 


(k) 


\-K{3]{k)  0 
- Kf\k )  1 

-^3(3)(^)  0 


0 

0 

1 


v3(k,k- 1), 


where  /3  is  the  3x3  identity  matrix,  and  v3  (k,k  —  \  )  was  defined  above. 

2. 6. 8  Track  Coasting 

As  discussed  in  Section  2.5.8,  when  a  track  is  to  be  coasted,  the  filter  gain  is  set  to 


(123) 


(124) 


K,  (k) 


0 

0 

0 


(125) 


With  the  Kalman  gain  set  to  zero,  the  remaining  equations  for  filtered  estimate  and  variance  are  the 
same  as  those  defined  above.  In  this  case,  the  filtered  estimate  is  given  by  the  prediction 


x3(k)  =  x3{k,k- 1) 


(126) 


or,  equivalently, 
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b3(k)  =  b3(k,k  -1) 

<  b3(k)  =  b3(k,k -1)  (127) 

b3(k)  =  b3(k,k- 1), 

and  the  filtered  variance  takes  the  form 

v3(k)  =  v3(k,k  - 1).  (128) 


3.  THE  ELEVATION  FILTER 

The  following  describes  the  mathematics  for  the  recursive  elevation  filter.  Both  the  initialization  and 
the  update  of  the  elevation  filter  are  described.  The  purpose  of  the  filter  is  to  estimate  elevation  and  to 
compute  time- varying  elevation  gates. 

The  elevation  filter  is  a  second-order  Kalman  filter  that  models  both  elevation  and  elevation  rate.  This 
same  filter  structure  is  used  for  the  second-order  bearing  filter.  The  filter  assumes  that  the  measurement 
error  variance  is  properly  defined  for  elevation  angle  as  a  function  of  measured  frequency  and  SNR. 


3.1  Overview 

The  elevation  filter  is  a  second-order  Kalman  filter.  The  filter  equation  corresponds  to  an  alpha-beta 
filter  and  is  given  by 


x(A:)  =  (129) 


where  for  the  elevation  filter,  we  define 


b(k) 

Kk) 


(130) 


0(k,k- 1) 


1  At 


(131) 


and 


z{k)  =  em(k). 


(132) 
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h  (*)=[i  o]. 


(133) 


In  the  above,  x(A)is  a  2  x  1  vector  of  filtered  estimates  valid  at  time  k- 1  from  the  previous  iteration 

and  valid  before  the  current  elevation  data  is  received.  At  is  the  difference  between  the  current  cluster 
measurement  time  and  the  last  time  the  track  was  updated  (time  of  previous  cluster  update). 
z{k)  =  em  (k)  is  the  elevation  angle  measurement  from  the  current  CD.  O  ( A ,  A  —  I )  is  a  2  x  2  state  or 

plant  transition  model  and  describes  the  emitter  motion  from  time  segment  A-l  to  k  .  K  (A)  is  a  2  x  1 
vector  of  Kalman  gain  coefficients  and  //(A)  is  a  1  x  2  measurement  vector. 

The  above  equation  for  the  filter  update  may  also  be  written  as 


x(A)  =  x(A,A-l)  +  //(A)  z(A)-//(A)x(A,A-l)]  . 


(134) 


where  the  above  is  written  in  terms  of  the  predicted  estimate 


x(k,  k- 1)  =  e(k,  k  - 1) . 


(135) 


Expanding  the  above  yields  the  two  equations  for  elevation  and  elevation  rate  as 


(A)  =  e(A,A-l)  +  //(1)  (A)[z(A)-z(A,A-l)] 
(A)  =  e(A,A-l)  +  //(2)  (A)[z(A)-z(A,A-l)], 


(136) 


where  (k)  is  equivalent  to  the  gain,  a,  used  by  an  alpha-beta  filter  and  /<" Ul  ( A )  is  equivalent  to 


(2)( 


(3/At. 


In  addition,  the  predicted  estimates  presented  above  are  defined  as  follows: 


e(A,A-l)  =  e(A-l)  +  Af  •e(A-l) 


(137) 


is  the  predicted  elevation  angle, 


e(A,A-l)  =  e{k- 1) 


(138) 


is  the  predicted  elevation  angle  rate,  and 


z(k,  k- 1)  =  H(k)  0>  ( k ,  k  - 1  )x(k  - 1)  =  e(k,  k  - 1) 


(139) 
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is  the  predicted  elevation  cluster  measurement. 

The  sequence  z(k)-z(k,k- 1)  is  known  as  the  innovation  sequence  and  represents  the  new 
information  to  the  filter. 

3.2  Initialization 

The  elevation  filter  may  be  initialized  using  one  of  three  equivalent  methods. 

Method  1 

In  this  method,  the  bearing  filter  covariance  matrix  is  initialized  with  large  values  on  the  diagonal: 


where 


-(0)  = 


and  ve  (0)  =  1012  I2 , 


(140) 


v,(0)  is  the  initial  2x2  filtered  covariance  value  and  fi  is  a  2  x  2  identity  matrix.  A  large  value  simply 
indicates  no  a  priori  information. 

Method  2 

In  this  method,  the  elevation  filter  is  initialized  from  the  first  cluster  measurement  of  elevation  angle. 
In  this  case,  we  define 


eJX) 

0 


and  ve(l) 


U1)  0 

n  i  niz 


(141) 


where  em(X)  is  the  elevation  measurement  provided  with  the  first  cluster  and  vce  ( I )  is  the  measurement 

accuracy  associated  with  the  elevation  measurement.  The  filter  iteration  starts  at  time  k  =  2  with  the 
second  measurement. 

Method  3 

In  this  method,  the  elevation  filter  is  initialized  from  the  first  two  cluster  measurements  of  elevation 
angle.  In  this  case,  we  define 
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em(2) 

(2)  ~  (!) 

At 


and  ve(2) 


Vce(2) 

Vce(2) 

At 


Vce(2) 

At 

2v„(2) 

(A,)2 


(142) 


The  filter  iteration  starts  at  time  k  =  3 .  This  is  the  preferred  method  since  it  starts  the  filter  when  a  CD 
elevation  measurement  has  correlated  with  a  first  point.  Therefore,  the  initial  elevation  CD  is  not  the  first 
ETF  track  entry. 

3.3  Predicted  Estimate 

When  a  new  CD  is  associated  with  the  elevation  track,  the  first  step  in  the  filter  update  process 
requires  that  the  ETF  estimate  be  predicted  to  the  time  of  the  current  associated  CD  (time  alignment).  For 
the  elevation  filter,  the  predicted  estimate  is  obtained  from  the  transition  matrix  defined  above,  the  filtered 
estimate  x(k ) ,  and  the  time  difference  between  the  associated  CD  and  the  ETF  time  (time  of  the  last 
track  update)  as 


x(k,k-\)  -  ®(£, &-1)jc(£-1) 


e[k,k- 1) 

e(£-l)  +  At-e(k-X) 

1 

V* 

i 

i _ 

1 

1 

_ 1 

(143) 


where  x{k,k- 1)  is  the  prediction  of  x(k  - 1 )  to  the  current  CD  time  and  x(k  —  l)  is  the  filtered 
estimate  obtained  from  the  ETF. 

3.4  Predicted  Variance  and  Plant  Noise  Estimates 


Once  the  prediction  is  computed,  computation  of  the  Kalman  filter  gain  is  required  to  merge  the 
elevation  CD  data  with  the  elevation  track.  The  Kalman  gain  computation  requires  an  estimate  of  the 
predicted  variance  and  the  measurement  variance.  The  predicted  variance  is  given  by 

ve(k,k-  l)  =  ve(*-l)  +  ve(*),  (144) 

where  ve  (k )  is  the  plant  noise  estimate.  For  the  elevation  filter,  the  plant  noise  is  set  to 


Ve(k)=be(k)-We-be{k)'  ’ 


(145) 


where 


bM) 


At 


2 

At 


and  w  =  <7 1 . 


(146) 
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a2  is  a  constant  and  should  be  approximately  set  to 

cre2  =0.1-^-/180/40  =  0.00005. 


Therefore, 


,  (k)  =  0.00005- 


f  At2  ^ 


V  2  J 

At 3 

~Y 


At 

2 

At2 


(147) 


(148) 


3.5  Filter  Gain  and  Measurement  Noise  Estimate 

The  Kalman  gain  is  given  by 


*(*) 


v,(k,k-\)H' 

Hv,(k,k-\)H'*vJkY 


(149) 


We  note  here  that  a  simple  divide  is  used  and  that  no  matrix  inverse  operation  is  required  since  the 
term  Hve(k,k -\)H'  is  a  scalar.  However,  to  compute  this  gain,  an  estimate  of  CD  elevation 

measurement  noise  is  required.  If  the  CD  elevation  measurement  noise  is  not  included  in  the  cluster 
message,  it  can  be  computed  from  the  single  pulse  specification  variance  according  to  the  equation 


v„(*)  = 


Vme(k) 


(150) 


where  vce  ( k )  is  the  measurement  noise  variance  in  elevation  associated  with  the  cluster  estimate, 
Vme{k)  is  the  single  pulse  measurement  accuracy  (equivalent  to  the  PDW  accuracy)  of  the  elevation 
angle,  and  n  is  the  number  of  pulses  reported  in  the  cluster  descriptor  message.  The  computation  of 
vme  ( k )  was  discussed  previously. 


The  Kalman  gain  may  also  be  written  in  terms  of  the  innovation  sequence  as 
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K(k) 


v,(k,k- l)-H' 
S(k) 


(151) 


where  S  (  k  j  is  the  scalar  covariance  of  the  innovation  sequence  defined  above  and  is  given  by  the  scalar 


S(k)  =  Hv,(k,k-l)H'  +  vJk). 


(152) 


3.6  The  Filtered  Estimate 

After  the  Kalman  gains  are  computed,  the  new  filtered  or  ETF  track  estimate  is  given  by 

x(k)  =  x(k,k- 1)  +  K(£)[z(£)  —  z{k,k  -1)] ,  (153) 

where  x(&)  is  the  filtered  elevation  estimate  and  z  ( k )  is  the  current  elevation  cluster  measurement.  As 
mentioned  above,  this  is  equivalent  to  computing 


e(£)  =  e[k,k  -l)  +  A(1)  (&) [z (&)-£(&,  A: -l)] 
e(£)  =  e(k,k-\}  +  K^'1  (A:)[z(A:)-z(A:,A:-l)]. 


(154) 


3.7  Filtered  Variance 

The  filtered  variance  of  this  estimate  is  simply 

ve(k)  =  \_I2-K(k)H~\ve(k,k-i),  (155) 

where  /0  is  the  2x2  identity  matrix,  and  ve[k,k  —  1 )  was  defined  above. 

3.8  Track  Coasting 

To  coast  a  track,  the  filter  gain  is  set  to 

0 
0 


(156) 
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With  the  Kalman  gain  set  to  zero,  the  remaining  equations  for  filtered  estimate  and  variance  are  the  same 
as  those  defined  above.  In  this  case,  the  filtered  estimate  is  given  by  the  prediction 


x(k)  =  x(k,k- 1)  (157) 


and  the  filtered  variance  takes  the  form 


ve(k)  =  ve(k,k-l) .  (158) 


4,  THE  ESM  PARAMETERS  FILTER 

The  ESM  parameters  filter  is  a  simple  first-order  Kalman  filter  that  may  be  used  to  update  both  agile 
and  non-agile  parameters  as  well  as  automatically  generate  gates  for  ESM  front-end  sorting.  The  filter 
equation  for  all  ESM  parameters  is  identical  to  an  alpha-only  filter  and  is  given  by 


x[k)  =  x(k-\)  +  a  z(k)-x(k- l)J, 


(159) 


where  x  is  the  parameter  to  be  stored  in  the  track  file  and  may  be  represented  by  frequency,  pulsewidth, 
pulse  repetition  interval  (PRI),  or  amplitude,  and  z  is  the  current  measured  and  averaged  parameters  from 
the  cluster  descriptor. 


4.1  Initialization 

Each  ESM  parameters  filter  may  be  initialized  using  one  of  two  equivalent  methods. 


Method  1 

The  filter  is  initialized  with 


x(0)  =  0  and  v/.(0)  =  1012 
k  =  0 

Gate  =  ESM  parameter  default  size, 

where  v  f  ( k )  is  the  associated  filtered  covariance  value.  A  large  value  simply  indicates  no  a  priori 
information.  In  this  case,  vf  (k)  is  a  scalar  since  the  filters  are  one-dimensional.  Gate  is  the  default  ESM 

gate  size  associated  with  the  front-end  default  gates  for  the  parameter  of  interest.  The  default  gate  results 
in  initial  clusters,  which  represent  first  points  to  the  ETF. 
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Method  2 

In  this  method,  the  ESM  parameters  filter  is  initialized  from  the  first  cluster  measurement.  In  this 
case,  we  define 


x(l)  =  zm(l)  and  v/(1)  =  vcz(l),  (160) 

where  zm(  1)  is  the  ESM  parameter  measurement  of  interest  (frequency,  PRI,  pulsewidth,  or  amplitude) 
and  vcz  (1)  is  its  associated  scalar  covariance.  vc,  ( k )  is  the  measurement  accuracy  associated  with  the 
ESM  parameter  measurement  provided  in  the  cluster  descriptor.  The  filter  iteration  starts  at  time  k  =  2  . 

4.2  Predicted  Estimate 

For  the  ESM  parameters  filter,  the  predicted  estimate  is  simply  the  filtered  estimate  since  no  higher- 
order  derivatives  are  used  for  ESM  parameter  tracking: 

x[k,k-\)  =  ;t(£-l),  (161) 

where  x[k,k- 1)  is  the  one  step  prediction  of  x(k- 1)  . 

4.3  Predicted  Variance,  Measurement  Noise,  and  Plant  Noise  Estimates 

The  predicted  variance  is  given  by 

Vp(k,k-\)  =  vf(k-\)  +  vp(k),  (162) 

where  v  (k  )  is  the  plant  noise  estimate.  For  the  case  with  jitter,  the  cluster  plant  noise  may  be  estimated 
as 


v«= 


(x  —x 

\  max  min  , 

n 


(163) 


where  v  (  k  )  refers  to  the  measured  plant  noise  or  spread  (max  minus  min)  obtained  from  the  parameter 

measurements  provided  in  the  cluster  descriptor,  and  (xmax,xmin)  are  the  largest  and  the  smallest 

parameter  measurement  in  the  cluster  descriptor,  respectively.  The  average  plant  noise  may  be  obtained 
by  using  an  appropriate  alpha  filter  with  gain  limited  to  0.1,  or 
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v,W  =  v,(*-l)  +  0.1. 


(  Xmax  Tulin  ) 

12 


-*,(*-  1) 


(164) 


where  this  filter  is  initialized  with  the  first  measurement  of  spread  provided  in  the  cluster  descriptor. 

This  will  accurately  estimate  the  plant  noise  if  the  measurement  noise  is  negligible  by  comparison 
(always  the  case  for  a  jittered  measurement).  If  the  estimate  of  parameter  spread  is  small,  then  no  agility 
is  present  and  the  plant  noise  should  be  set  to  one-tenth  of  the  measurement  noise  for  heavier  smoothing. 
This  is  checked  for  by  implementing  the  test  described  below. 


For 


vpc  (k)  <  4vc.  (£)  or  Freq  Type  = 'Fixed' , 
set 


10 


4.4  Filter  Gain 

The  filter  gain,  alpha,  is  a  function  of  the  predicted  variance.  This  gain  is  computed  after  the  cluster  is 
received  and  begins  the  filter  iteration  for  the  current  update.  The  Kalman  gain  is  given  by 


a(jA-  Vp(k’k~l) 

1  ~  vp{k,k-\)  +  vcz{k)' 


(165) 


However,  to  compute  this  gain,  an  estimate  of  parameter  measurement  noise  is  required.  This  is 
computed  as 


v„(*)  = 


vmAk) 


(166) 


where  v  „  (£)  is  the  measurement  noise  variance  associated  with  the  cluster  estimate,  vm_  ( k )  is  the 

single  pulse  (PWD)  measurement  accuracy  of  the  parameter  of  interest,  and  n  is  the  number  of  pulses 
reported  in  the  cluster  descriptor  message. 

If  unknown,  the  cluster  measurement  variance  may  be  approximated  by 


v 


cz 


(1) 


(1.5  -res) 
n 


(167) 


where  res  is  the  resolution  of  the  parameter  of  interest. 
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4.5  Filtered  Estimate 

The  filtered  estimate  is  given  by 

x(A')  =  x(A:-l)  +  a(^)[z(^)-x(A:-l)],  (168) 

where  x  ( k  j  is  the  filtered  ESM  parameters  estimate  and  z  ( k  j  is  the  current  ESM  cluster  measurement. 
Due  to  the  initialization,  for  the  first  ESM  parameter  update  this  should  yield 


a  =  1  and  x(l)  =  z(l) .  (169) 


4.6  Filtered  Variance 

The  filtered  variance  of  this  estimate  is  simply  given  by 


vf(k)  =  [l-a(k)~\vp(k,k-l).  (170) 

For  the  first  cluster  update,  this  will  be  approximately  equal  to  the  initial  cluster  measurement  variance 


vf  (l)  =  vA1)- 


(171) 


5.  SUMMARY 

An  ESM  angle  tracker  was  developed  to  provide  highly  accurate  angle  and  angle  rate  estimates  to  a 
combat  system.  The  tracker  was  developed  to  operate  with  the  NRL  ESM-ATD  that  provides  accurate  but 
intermittent  bearing  and  elevation  reports.  This  report  describes  the  theory  behind  the  bearing,  elevation, 
and  ESM  parameters  filters  contained  within  the  tracker. 

The  bearing  filter  is  an  IMM  filter  that  automatically  adjusts  to  optimally  track  stationary,  constant- 
velocity,  or  constant-acceleration  targets.  The  IMM  bearing  filter  uses  three  standard  Kalman  filters  of 
different  orders  as  a  basis  for  filtering  and  predicting  emitter  bearing  angle.  Details  on  how  to  combine 
estimates  from  the  first-,  second-,  and  third-order  Kalman  filters  within  the  IMM  are  developed  and 
described  in  this  report. 

The  inclusion  of  ESM  parameters  in  the  tracking  algorithm  makes  it  possible  to  associate  and 
correlate  intermittent  bearing/elevation  measurements  with  existing  tracks  in  dense  scenarios.  The  ESM 
parameters  filters  are  first-order  Kalman  filters,  but  they  have  the  ability  to  track  both  constant  and 
jittered  parameters  by  automatically  adjusting  the  filter  plant  noise.  The  technique  for  estimating  plant 
noise  from  the  received  measurements  is  described  in  this  report. 
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